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We show that topological vortex pumping can be implemented for a dilute Bose-Einstein conden- 
sate confined in a magnetic time-averaged orbiting potential trap with axial optical confinement. 
Contrary to earlier proposals for the vortex pump, we do not employ an additional optical potential 
to trap the condensate in the radial direction, but instead, this confinement is provided by the 
magnetic field throughout the pumping cycle. By performing numerical simulations based on the 
spin-1 Gross-Pitaevskii equation, we find that several pumping cycles can be carried out to produce 
a highly charged vortex before a majority of the particles escapes from the trap or before the vortex 
splits into singly quantized vortices. On the other hand, we observe that an additional, relatively 
weak optical plug potential is efficient in preventing splitting and reducing particle loss. With this 
study, we aim to bring the vortex pump closer to experimental realization. 
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I. INTRODUCTION 

Bose-Einstein condensation in alkali gases was ob- 
served experimentally in 1995 [H-Q. A few years later, 
these pioneering experiments were followed by the cre- 
ation of singly quantized vortices 0, Q and vortex lat- 
tices in such systems. Since then, the study of 
vortices in Bose-Einstein condensates (BECs) has flour- 
ished both theoretically and experimentally [T3, E| due 
to their close connection with phase coherence and super- 
fluidity. Especially, their stability has been the subject 
of extensive research fl2l - |24| . 

In a loop encircling a quantized vortex, the phase of 
the condensate order parameter undergoes an integer 
multiple k of 2tt windings. In principle, a vortex in a 
BEC can have any winding number k. However, it is 
well known that a vortex with |«| > 1 typically has a 
higher energy than the corresponding number of sep- 
arated single-quantum vortices. Consequently, vortices 
with large winding numbers are prone to splitting pBT - 
|32| . which renders them challenging to create with dy- 
namical methods, such as using a focused laser beam to 
stir Q or slice through [33| the BEC, rotating it with an 
asymmetric trap potential [Hj], or colliding condensates 
separated by tailored optical potentials j35(. Being able 
to produce vortices with large winding numbers would 
provide access to novel vortex splitting patterns bey ond 
the typical linear chain that prevails for < 4 pa l30j. 
Due to the distinct nature of the different splitting pat- 
terns predicted for large values of k [HI, observing the 
decay of such vortices would allow for a lucid compar- 
ison between theory and experiment. Moreover, it has 
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been speculated that giant-vortex splitting may create 
necessary conditions for the initialization of superfluid 
turbulence [H, |37[. 

In addition to the above-mentioned dynamical meth- 
ods, vortices can be created with the so-called topological 
phase engineering technique (38l-[4l| (see Ref. [42| for re- 
view), in which the spin degree of freedom of the BEC is 
controlled adiabatically by a time-dependent nonuniform 
magnetic field. Since the method does not rely on the 
relaxation of condensate dynamics, it is especially well 
suited for producing multiquantum vortices. Indeed, the 
first two-quantum and four-quantum vortices in dilute 
BECs were created by applying the technique to spin-1 
and spin-2 BECs confined in a magnetic Ioffc-Pritchard 
(IP) trap E3. 

Subsequent theoretical studies have demonstrated that 
the topological phase engineering technique can also be 
used to implement a so-called vortex pump j44TJ48j . In 
this device, a fixed amount of vorticity is added to the 
BEC in each control cycle, and thus its repeated applica- 
tion would — stability issues notwithstanding — enable the 
creation of vortices with arbitrarily large winding num- 
bers. The original proposal |44j involved creating 2F 
quanta of vorticity per cycle in a spin-i* 1 BEC with a 
magnetic field configuration consisting of the standard IP 
trap and an additional hexapolc magnetic field. Backed 
by numerical simulations, the pump was shown to be 
operable both fully adiabatically and partly nonadiabat- 
ically. Later, Xu et al. (45| presented a different pumping 
cycle for the IP trap in which the hexapole field was re- 
placed with a uniform transverse field. Unfortunately, 
both of these control cycles suffer from the fact that the 
magnetic fields provide radial confinement only during 
part of the cycle, and thus, the fully adiabatic operation 
of the pump necessitates an optical trap to confine the 
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BEC radially. Since the purpose of the IP trap has been 
to confine the atomic cloud in the first place, there has 
been little incentive to augment it with an optical trap 
or a hexapole field. Therefore, pumping schemes not re- 
quiring such extra ingredients would be desirable from a 
practical standpoint. 

Recently, Xu et al. showed that vortex pumping can be 
applied to quantum superpositions to generate counter- 
circulation states [I?} . The authors used a novel control 
cycle which is particularly suitable for the time-averaged 
orbiting potential (TOP) trap and has the advantage 
that the radially confining quadrupole field can be kept 
on throughout the entire cycle. However, since only one 
of the components in the superposition state could be 
trapped magnetically, a three-dimensional optical trap, 
as well as a strong optical plug potential piercing the 
vortex core, had to be employed [471 ] . 

The aim of this article is to bring the vortex pump 
closer to experimental realization by showing that it can 
be implemented with mature, existing technologies al- 
ternative to the ones considered in Refs. [4J-[48|]. To 
this end, we demonstrate that vortices can be efficiently 
pum ped in the TOP trap without using, in contrast to 
Ref. [42j , additional optical potentials to confine the BEC 
in the radial direction or to pin the vortex core. In- 
stead, the radial confinement is provided solely by the 
magnetic field throughout the pumping process, and op- 
tical trapping is required only in the axial direction. We 
present simulations based on the Gross-Pitaevskii equa- 
tion which indicate that several pumping cycles can be 
carried out before a majority of the particles escapes from 
the trap or before the generated multiquantum vortex 
splits into sin gly quantized vortices due to dynamical in- 
stabilities [32u49j . On the other hand, we also show that 
even a relatively weak optical plug potential is efficient in 
preventing the splitting and in reducing the loss of parti- 
cles, thereby enabling the controlled creation of isolated 
vortices with large winding numbers. 

The remainder of this article is organized as follows. 
In Sec. |nl we present the zero-temperature mean-field 
theory of the spin-1 BEC, describe the control cycle of 
the vortex pump, and discuss in detail the confinement 
of the condensate during the cycle. Section Hill presents 
our numerical results, which we relate to realistic exper- 
imental setups in Sec. IIVI Finally, Sec. [V] concludes the 
article with a discussion. 



II. THEORY AND METHODS 

A. Mean-field model 

We consider a dilute spin-1 BEC in the zero- 
temperature limit, thereby neglecting the possible effects 
due to noncondensed atoms. In the standard mean- 
field treatment, the spin-1 condensate is described by 
a three-component order-parameter field that we write 
in the cigenbasis of the spin-1 matrix F z as ^ = 



(^+1, ^o, ^-i). Its time dependence is given by the spin- 
1 Gross-Pitaevskii (GP) equation [HH l52j 

ihd t $(r,t) = (u + g n &$ 

+g s &F$-F)$(r,t). (1) 
The single-particle Hamiltonian operator H. is given by 

h 2 

H = -—V 2 + V opt (r)+ f i B g F B{t)-F, (2) 
Zm 

where m denotes the atomic mass, g F is the Lande fac- 
tor, fiB is the Bohr magneton, B(r, t) denotes the ex- 
ternal magnetic field, and F = (F x , F„, F z ) is a vec- 
tor of the standard spin-1 matrices |53f . Optical poten- 
tial terms are contained in V op t(r) = V tY (z) + V p \ ug (r), 
where V tl -(z) = muj 2 z 2 /2 is a strong axial harmonic 
trap and V^i ug (r) = Acxp (— r 2 /d 2 ) describes a possi- 
bly present Gaussian-shaped repulsive plug potential of 
amplitude A > and width d. Here, r = \/ x 2 + y 2 is 
the radial coordinate. The coupling constants g n and 
g s appearing in Eq. ([1} measure the strengths of the 
local density-density and spin-spin interactions, respec- 
tively. They are related to the s-wave scattering lengths 
and a*- 2 ) into spin channels with total spin and 
2h by the expressions <? n = Airh 2 (a^ 1 + 2et( 2 )) /3m and 
g s = inh 2 (a^ 2 ) — a*- -*) /3m. The order parameter is nor- 
malized such that J d 3 r\Ir ^ = No, where Nq is the num- 
ber of particles in the BEC. 

B. Magnetic fields and the pumping cycle 

The operation principle of the vortex pump is to con- 
trol the spin degree of freedom of the condensate locally 
by slowly tuning the magnetic field B(r, t) in a cyclic 
manner such that the system acquires a fixed amount of 
vorticity per cycle [4414481 ]. In the pumping scheme con- 
sidered here, the spin-1 atoms are assumed to be magnet- 
ically confined in the standard TOP trap . It consists 
of a quadrupole field B q , which has axial symmetry about 
the z direction, and a rapidly rotating, spatially uniform 
magnetic field B ro t oriented along the xy plane. In addi- 
tion, we assume that the TOP trap is augmented with a 
uniform axial bias field B^(t)z that can be controlled in- 
dependently of the other fields. The total magnetic field 
can be written as 

B(r, t) = B q (r) + B rot (t) + B 7 b (i)z, (3) 

where B q (r) = B' (xx + yy — 2zi) is the quadrupole 
field with the radial gradient B' and the rotating trans- 
verse bias field is given by 

B rot (*) = -B rot (t) [cos (cj rot t) ic + sin (w rot t) y] , (4) 

where w ro t denotes its angular frequency of rotation 
about the z axis. The bias field strengths are assumed 
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to be bound by B such that B^(t) 6 [— Bq,Bq] and 
B m t{t) G [0, So]. We point out that the field config- 
uration of Eg. (E0) has already been employed in BEC 
experiments [5J, [55| . 

In order to facilitate vortex pumping in the TOP trap, 
we use the control cycle presented in Fig. [1] It is car- 
ried out by tuning two parameters, B\ and B Iot , and 
can be divided into part A (0 < t < T A ) and part B 
{Ta < t < T A + T B =: T). Part A is similar to the origi- 
nal proposals of topol ogical p hase engineering [38U4ll Hj| 
and the experiments pa l43l. [56l - [59j . and it is responsible 
for increasing the circulation in the spin-1 BEC by two 
quanta. It is executed by reversing the axial bias field 
with the rotating field switched off, 



(a) 



Brot(i)=0, 

B*(t) = B'potan 



2t-Ti 



Ti 



arctan 



\B'p J 



< i < T A , 
(5) 

where Bq should be large enough to render the BEC es- 
sentially spin polarized along the z axis at t = and 
t = Ta- To improve adiabaticity, the time dependence 
for B\ has been chosen such that spins at a distance of 
po from the z axis are turned with constant speed, but 
part A can also be performed by reversing B^(t) linearly 
in time (43|. In part B, the axial bias field is returned to 
its initial value while ramping up and down the rotating 
field, 



' B Iot (t) = B sw/3(t), 
^B h z {t) =B o cos0(t), 



T A <t<T, 



(0) 



where (3{t) = it (t — Ta) /Tb. Part B was originally pro- 
posed by Xu et al. [47| , and it is designed to preserve the 
accumulated vorticity. The cyclic repetition of parts A 
and B will therefore increase the vortex winding number 
of the spin-1 BEC by two per cycle. 

To efficiently steer the condensate spin by the magnetic 
field B (r, t) requires that the Zeeman energy dominates 
over the kinetic energy at each point in space. Hence, 
to guarantee adiabaticity, |B (r, t) | should be sufficiently 
large in the region occupied by the BEC. This condi- 
tion is not fulfilled at the origin when B\ crosses zero 
at t = TX/2, and thus it is desirable to prevent parti- 
cles from entering this area. This can be accomplished 
by introducing the repulsive plug potential V^i ug (r) along 
the z axis. The plug not only improves adiabaticity but 
also serves to stabilize the created multiquantum vortex 
against splitting [H, H^] . The plug can be realized for 
pancake-shaped BECs by a focused blue-detuned laser 
beam as has been done in various experiments [3l.l8ll9l.l60l- 
|62| . In this article, we present results for vortex pumping 
both with and without the plug potential. 



C. Confinement during pumping 

An essential difference between this article and ear- 
lier work concerning the vortex pump (43 - [48| is that 



o 

^0 
cq 





FIG. 1. (a) Control cycle of the vortex pump in the (B r ot, ) 
plane, where B Iot and £>!? are the strengths of the transverse 
and axial magnetic bias fields, respectively. The cycle starts 
at (Brot,Bz) = (0, — -Bo) and proceeds clockwise, (b) Time 
dependence of B rot (solid line) and B\ (dashed line) [Eqs. ([SJ) 
and ([6])] during the cycle. The cycle is divided into parts A 
and B as indicated. In part A, the axial bias field B\ can also 
be reversed linearly. 



here we never employ an optical trapping potential in 
the radial direction. Instead, this confinement is pro- 
vided by the magnetic field throughout the entire pump- 
ing cycle. In the case of spin-I BECs, the magnetically 
trapped weak-field seeking state (WFSS) corresponds to 
the highest-energy eigenstate of the Zeeman Hamilto- 
nian c/j^bB ■ F, with the effective trap potential given 
by the spatially varying eigenvalue |<7i?^BB(r, t)\. Even 
though this potential is solely responsible for the ra- 
dial confinement, a strong optical trap is still needed 
in the axial direction to keep the atomic cloud centered 
around z = throughout the pumping cycle. Changes 
in £?k shift the z coordinate of the zero-value point of 
the total magnetic field, and without the optical z con- 
finement, carrying out the cycle would merely move the 
whole BEC along the z axis. Hence, the Hamiltonian 
in Eq. ^ includes the harmonic axial trapping potcn- 
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tial Vtr(z) = muj 2 z 2 /2 with a trap frequency uj z that 
is assumed to be large enough to render the condensate 
pancake shaped, which means that the order parameter 
can be taken to have the form (r, t) = "J^d (x, y, t) ((z), 
where ((z) = exp (— z 2 /2a 2 ) / -^/"faf and a z = y/h/mu) z 
is the axial oscillator length. This enables us to inte- 
grate out the z variable in Eq. ([1} and obtain an effec- 
tively two-dimensional GP equation with the magnetic 
field determined at z = 0. 

Let us consider the shape of the magnetic potential in 
the vicinity of the origin. During part A of the cycle 
(0 < t < Ta), the magnetic field strength is given by 
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(B'y) 2 

2B'\B h \ 
! £1 z 4- 



B' 2 
2\W\ 



2B'z] 



where in the expansion we have neglected third- and 
higher-order terms in B'r/\B^\ and There- 
fore, the magnetic field at t — gives rise to an approx- 
imately harmonic potential in the radial direction with 
the trap frequency 



oj = B'(0) ] 



9fPb 



mB 2 b (0) 



(8) 



It is convenient to measure all quantities in terms of loq, 
and thus we express lengths in units of the corresponding 
oscillator length ao = \J h/mwo ~> a z , energies in units 
of Swo; time in units of 1/loq, and the magnetic field in 
units of ^ojo/|3f|mb- Variables expressed in these units 
arc henceforth denoted with a tilde. 

According to Eq. ([7]) , the profile of the radial confine- 
ment will change during part A. Initially, the potential is 
harmonic, with the effective trap frequency oc |S^(t)| -1 / 2 
increasing in time with decreasing |.S°(£)|. At t = Ta/2, 
the trap becomes purely linear in r with the gradient 
\gp\[iftB' . The axisymmctric modulations of the trap 
profile will cause shrinking of the BEC and excitation of 
its breathing mode. Although these effects do not criti- 
cally hinder the operation of the pump, they can never- 
theless be reduced by introducing time dependence into 
B' such that it is decreased when \B Z \ is ramped down 
during part A. Here, we use the dependence 



B'(t) = 



B' 
BL 



B' 



2t-T A 



,0<t<T A , 

T A < t < T. 
(9) 

with B' = B'(0) and B' min tn 0AB' o . Varying B' is not 
necessary but it improves the accuracy of the pump and 
reduces the loss of particles. 

In part B, the frequency of the rotating bias field, w ro t, 
is chosen to be low compared with the frequencies of tran- 
sitions between different magnetic substates but large 
compared with the effective radial trap frequency. Typi- 
cally, magnetic trap frequencies are of order 10 2 Hz while 
the transition frequencies are of order 10 6 Hz. Thus, a 
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reasonable choice would be, e.g., w ro t ~ 10 4 Hz. These 
conditions ensure that the atoms will not undergo tran- 
sitions to other substates and be lost from the magnetic 
trap but instead move in an effective potential given 
by the time average of the instantaneous magnetic po- 
tential over one rotation period of the field B rot . For 
TA<t<T, the fast-time-averaged field strength can be 
expanded in a power series in B't/Bq and B'\z\/Bq as 

^ / dt'\B{t')\ ^ B Q -2B'z cos (3 + ^- (l + cos 2 /3)r 
2tt Jt, 4±Sq 



2B' 2 sm 2 B , 

+ 1T^ Z ' 

-DO 



(10) 



where the third- and higher-order terms have been 
discarded and the interval of intervation is It = 
[t — ir/tdrot, t + 7r/cj ro t]. Equation (fTU| implies that the 
effective radial trap frequency will decrease by ~ 29% 
during part B, but this should not significantly disturb 
the pumping process. 



III. RESULTS 

We study the temporal evolution of a spin-1 BEC dur- 
ing vortex pumping by numerically solving the GP equa- 
tion, Eq. (|T|), with the time dependence of the mag- 
netic field B(r,f) given by Eqs. ©, ©, and ©. Af- 
ter factoring out the z dependence of the order pa- 
rameter as \&(r, i) = (x, y, t) ((z), Eq. ((T|) is dis- 
cretized on a uniform grid with a finite-difference method 
and integrated in time with a split-operator approach. 
The dimcnsionlcss coupling constants are chosen to have 
the values g n = N mg n / y / 2nh' l dl = 250 and g s = 
N mg s / y / 2Trh 4 a 2 = — 0.01t/ n , the latter correspond- 
ing to spin-1 condensates of 87 Rb [63 66]. The dura- 
tions for parts A and B of the control cycle are given 
by Ta = uiqTa = 3 and Tb = ojqTb = 2, respec- 
tively, and the parameters in Eqs. ([5]), ([5]), and © 
have the values Bq = \gF\fJ.BBo/hcuo = 200, po = 
po/a = 5, B' = \g F \p B a B' /hwo = y/200, and B' min = 
\gF\pBaoB' min /Huio = 6. The Lande factor g F is taken to 
be negative as in the case of spin-1 87 Rb. The frequency 
of the rotating field is set to (D rot = uj Iot /ujo = 85. We 
present results both with and without an optical plug 
potential of amplitude A = A/Hluq = 10 and width 
d = d/ao = 2. Before the pumping is started, a re- 
laxation method is used to bring the BEC to the lowest- 
energy WFSS with the magnetic field in its t = config- 
uration. 

Figure [2] shows the squared moduli and the complex 
phases of the most relevant order-parameter components 
*+i (at t = IT, I e N) and (at t = IT + T A ) during 
the pumping process. The accumulation of two quanta of 
vorticity during part A of each cycle is clearly visible in 
the phase fields at t = IT+Ta- Part B of the cycle, during 
which the rotating bias field is on, is observed to leave the 
vorticity unaffected. The pumping also causes breathing 
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FIG. 2. (Color online) Areal particle density and complex phase of the order-parameter components SP+i (at t = IT, I £ N) and 
9-i (at t = IT + Ta) in the x?/ plane during pumping (a) without and (b) with an optical plug potential (amplitude A = Whwo, 
width d — 2ao). The field of view in the panels is (a) 12ao x 12ao and (b) 14ao x 14ao. The two parts of the control cycle have 
the durations Ta = 3/wo and Tb = 2/ljo, and the dimensionless coupling constants are set to g n = 250 and g s = —0.01g n . 



of the BEC, as indicated by its oscillating spatial extent 
and by the nonzero radial derivatives of the phase fields. 
The excitation of the breathing mode is attributed to 
the changing magnetic confinement during the pumping 
cycle [see Eqs. © and ([TO])]. 



Axisymmetric vortex states with large winding num- 
bers k have been found to be dynamically unsta- 
ble against splitting in pancake-shaped, harmonically 
trapped single-component BECs for most values of the 
interatomic interaction strength, with the degree of in- 
stability generally increasing with increasing k (24l . |32| . 
Therefore, when the stabilizing plug potential is not em- 
ployed in the pumping, the created multiquantum vor- 
tex is expected to split after it has accumulated a suf- 
ficiently large winding number. In Fig. (2Ja), the onset 
of splitting is visible around t = T + Ta, when k = 4. 
As shown for t = 5T + Ta, the process eventually re- 
sults in a line of singly-quantized vortices and is thus 
ascribed to a dynamically unstable excitation mode with 



orbital angular momentum of ±2h per particle with re- 
spect to the condensate [3^|. On the other hand, when 
the relatively weak plug is employed [Fig. [3jb)] , the vor- 
tex does not split despite its significant breathing, and 
a nearly symmetric 12-quantum- vortex state is observed 
at t = 5T + Ta- We have confirmed numerically that 
the plug amplitude A can be subsequently ramped down 
without destroying the state. 

Due to the finite pumping period T and magnetic field 
strength |B|, the pumping process is not perfectly adi- 
abatic and there are spins that do not follow the lo- 
cal magnetic field. Since these spin components are no 
longer trapped, they escape the condensate region. Con- 
sequently, the number of particles in the trap decreases 
during the process. The loss rate depends on the degree 
of adiabaticity of the pump, i.e., on the pumping speed 
and on the local magnetic field strength |B(r)|. 

Figure [3] presents the number of particles in the trap, 
A" = J r<R d 3 r\5'\5 < Nq, and their average orbital angu- 
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lar momentum (L z ) /N = J r<R d 3 r& [z ■ (r x V)] * 
as functions of time during the pumping process. Here, 
the particles are considered lost after their distance from 
the z axis exceeds R = 13ao. Therefore, the curves in 
Fig. [3] also include contributions from unconfined atoms 
that have not drifted away from the trap region. Even 
after taking this into account, we find that a significant 
portion of the atoms still remains in the WFSS at t = 5T. 
As expected, the optical plug [Fig. E^b)] is observed to 
significantly reduce the loss of particles. 

Whereas the number of vortices in a scalar conden- 
sate is quantized, the orbital angular momentum (L z ) is 
a continuous quantity. The angular-momentum curves 
in Fig. [3] indicate the increment of vorticity by two dur- 
ing part A of each cycle and show that (L z ) /N increases 
monotonously during the first few cycles. The consider- 
able deviation of (L z ) from the ideal value 21HN after I 
cycles is mainly due to the contribution of the untrapped 
atoms that remain in the region r < R. Moreover, the 
slight increase in (L z ) /N during each part B is attributed 
to the small center-of-mass motion induced by the rotat- 
ing transverse bias field. 



IV. EXPERIMENTAL FEASIBILITY 

Let us briefly relate the proposed topological vortex 
pump to a realistic experimental setup. As an exam- 
ple of previously realized experimental parameters, the 
experiment of Ref. [62j studied condensates of 2 x 10 6 
87 Rb atoms in the 5 2 Sy 2 \F = 1, m F = -1) state. The 
atoms were confined in a TOP trap with a bias field of 
£? ro t = 5 G and a quadrupole field with the radial gradi- 
ent B' = 27 G/cm. An additional red-detuned Gaussian 
laser beam propagating in the xy plane provided strong 
confinement along the z axis and negligible confinement 
in the radial plane; the TOP trap provided the signifi- 
cant portion of radial confinement. This laser beam had 
a wavelength of 1090 nm, a power of about 0.5 W, and 
radii of about 20 /Ltm along z and 2 mm along r. The 
combined optical and magnetic trap gave trapping fre- 
quencies of (wo, w z) — 27r x (8, 90) Hz, low compared with 
the TOP trap rotation frequency of uj xo t = 2ir x 4 kHz. 

As a proposed implementation of the vortex pump, we 
consider the spin-1 condensates of 87 Rb and the following 
field parameters. First, in place of a red-detuned trap- 
ping laser, we assume the use of a blue-detuned beam 
that has a Gaussian profile in the xy plane and a first- 
order Hermite-Gauss profile along z. Since the atoms 
would be trapped in the dark region between the two 
halves of the beam, the laser would not provide any ra- 
dial confinement. A 1-W, 532-nm beam, with Gaussian 
radii of approximately 9 /im along z and 500 ^m along r, 
will give an axial trapping frequency of to z w 2tt x 900 Hz. 
This field remains constant throughout the pumping cy- 
cle and is large enough to support the atoms against grav- 
ity and the magnetic-field forces along z. 

Secondly, for the magnetic field parameters, we as- 



sume the field values given above that correspond to 
Ref. [62j . Thus, we take the maximum strength of the 
uniform fields to be B = 5 G [Eqs. ([5]) and ||6J], which 
is reached for the axial bias field B\ at the beginning 
and end of part A of the cycle, i.e., at t = and 
t = Ta- Between these times, the rotating component 
of the TOP trap is off, the quadrupole field with the gra- 
dient B' = 27 G/cm is on, and the bias field pushes the 
zero- value point of the total magnetic field B from nearly 
1 mm above the pancake-shaped BEC to nearly 1 mm be- 
low the BEC. This satisfies the approximation of Eq. ([7]) 
and the harmonic approximation. For this trap, Eq. ([5)1 
yields uiq rs 2tt x 11 Hz, and hence the simulations in 
Sec. IIIII correspond to Ta = 3/luq = 44 ms. As assumed 
in Eq. ([5]), B' (t) can optionally be ramped to a minimum 
value of B' min w 11 G/cm during part A. In part B, the 
rotating bias field B rot is ramped on while the bias field 
B\ reverses direction, bringing the zero- value point of B 
in a spiraling trajectory around the BEC prior to B lot 
being ramped back off. Our simulations fix the duration 
of this stage at Tb = 2/wo = 29 ms. 

Additionally, for the parameters assumed above, the 
implementation of an optical plug would be straightfor- 
ward. The assumed value of the plug radius used in our 
simulations, d = 2ao, corresponds to a Gaussian 1/e 2 
beam radius of 2y2ao ~ 9 /im, similar to the beam used 
in Ref. 

For numerical convenience, our simulations have as- 
sumed smaller numbers of atoms than would be ideally 
used in an experiment, as well as lower values of magnetic 
fields than those typically found in TOP traps that have 
trap frequencies of the order stated above [67| . Neverthe- 
less, based on the validity of the harmonic approximation 
and the readily achievable time and length scales, the pri- 
mary features seen in the simulations should be preserved 
and observable with experimentally feasible parameters. 



V. CONCLUSION 

In summary, we have discussed how to implement a 
vortex pump for a BEC in a TOP trap, resorting only 
to standard experimental techniques and magnetic field 
configurations that are already available in BEC labora- 
tories. We showed that the pumping can be carried out 
without using an additional optical potential to trap the 
atoms in the radial direction. Instead, the radial con- 
finement is provided solely by the magnetic field of the 
TOP trap, and a harmonic optical potential is employed 
only in the axial direction. Our simulations showed that 
even if the pumped multiquantum vortices are not stabi- 
lized by a Gaussian-shaped plug potential piercing their 
core, several pump cycles can still be carried out before 
the vortex splits clearly. On the other hand, already a 
relatively weak plug potential was found to prevent the 
splitting and to reduce the loss of atoms from the trap. 

The experimental realization of the vortex pump would 
represent an important milestone in vortex physics, since 




FIG. 3. (Color online) Particle number N and average orbital angular momentum {L z )/N of the BEC as functions of time 
during the vortex pump simulations presented in Fig. [2] (a) no optical plug potential; (b) Gaussian-shaped plug of amplitude 
A = lOhuo and width d = 2a r . 



it would provide a controlled method to produce almost 
any desired amount of vorticity. From a theoretical point 
of view, the vortex pump is a fascinating example of adi- 
abatic quantum dynamics for which the control param- 
eters of the system are varied cyclically but the system 
does not return to its initial eigcnspacc. In fact, the 
appearance of vortices can be interpreted as the accumu- 
lation of a position-dependent geometric Berr y ph ase (68j 
for individual spins of the condensate atoms [38|-|4J| . 
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